Detecting Non-Abelian Anyons by Charging Spectroscopy 
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Observation of non-Abelian statistics for the e/4 quasiparticles in the ~ ^ fractional quantum Hall state 
remains an outstanding experimental problem. The non-Abelian statistics are linked to the presence of additional 
low energy states in a system with localised quasiparticles, and hence an additional low-temperature entropy. 
Recent experiments, which detect changes in the number of quasiparticles trapped in a local potential well as 
a function of an applied gate voltage, Vg, provide a possibility for measuring this entropy, if carried out over 
a suitable range of temperatures, T. We present a microscopic model for quasiparticles in a potential well 
and study the effects of non-Abelian statistics on the charge stability diagram in the Va — T plane, including 
broadening at finite temperature. We predict a measurable slope for the first quasiparticle charging line, and an 
even-odd effect in the diagram, which is a signature of non-Abelian statistics. 

PACS numbers: 73.43.Cd, 05.30.Pr, 71.10.Pm 



The unambiguous observation of particles obeying non- 
Abelian statistics remains an outstanding experimental chal- 
lenge in condensed matter physics. The Moore-Read frac- 
tional quantum Hall state (FQH) believed to be realized 
at filling fraction u ^ 5/2, is one of the most promising can- 
didate phases to exhibit such quasiparticles (QPs) The 
Moore-Read state is predicted to support QPs of charge ±e/4; 
for N such QPs, localized and well separated from each other, 
there should be a nearly degenerate set of ground states, with 
multiplicity 2^/^^^. For temperature T larger than the split- 
ting of these ground states, but smaller than the gap to higher 
excited states, this degeneracy contributes an effective entropy 
to the system, the non-Abelian entropy. 

Non-Abelian statistics predicts that pairs of QPs can inter- 
act to form two distinct states, or fusion channels, f, com- 
monly denoted as / = 1, ip. Ina finite system, the two states 
have different energies, and the ground state is unique; for T 
below the splitting between the two, the non-Abelian entropy 
is lost. 




FIG. 1. Cartoon charge stability diagram, showing only peak centres 
(no broadening). Vertical axis is temperature T; horizontal axis /i is 
the chemical potential for charged QPs, controlled in experiment by a 
gate potential. Red dashed lines are for Abelian particles. Blue solid 
lines correspond to non-Abelian QPs in a tightly confining well. Ajv 
is the gap to excited states for particles, and sets scale for other 
entropic effects. Notice even-odd effect for non-Abelian anyons. 

There are several recent theoretical proposals for tech- 



niques to observe this entropy through bulk measurement 
of thermodynamic and transport properties |3-6|. Recent 
measurements in this direction of thermoelectric response at 
h' ~ 5/2 are encouraging |7|. These theoretical proposals 
assume that all QPs are well separated, such that degeneracy- 
lifting interactions are weak or non-existent. However, recent 
local electronic charge-sensing measurements, using a single- 
electron transistor (SET) |8|, suggest that QPs tend to trap in 
local potential wells due to electrostatic disorder, which may 
be tightly confining and contain more than one QP. Confined 
QPs split their degeneracy through two means: Majorana ex- 
change 19] Ho), present even for stationary QPs, and as we 
show here, an orbital splitting from interchange of the charged 
QPs, which can dominate in special cases. 

In this letter, we study the charging spectra of local quasi- 
particle traps. Such traps may be induced by disorder or de- 
fined by gates. Their spectra reflect the QP statistics, just as 
electronic dot spectra reflect the spin and fermionic statistics 
of electrons. We show that low-frequency SET charge-sensing 
measurements, which provide only thermally-averaged infor- 
mation regarding the dot spectra, are sufficient for extracting 
non-Abelian signatures. At low but experimentally accessible 
T, we predict a robust temperature evolution of the = — 1 
transition, and an even-odd effect in the evolution of the charg- 
ing spectrum for several non-Abelian anyons. This effect 
should be visible for T below the relevant gaps to excited 
states for A^ particles, which we calculate for A^ = 1,2. 

The experiments of [iS] measure the change in potential at 
the SET induced by a change 6Vg in the potential applied to 
a backgate on the sample. If there is a single disorder-induced 
well close to the SET, the measured signal is inversely propor- 
tional to the compressibility of the well, k = ^1,^^ , where fi 
is the QP chemical potential in the vicinity of the well. For 
an isolated well, the relation between 6Vg and the change 
in /i should be linear, but the constant of proportionality is 
geometry-dependent, as screening depends on the local en- 
vironment as well as the distance to the gate ifTTI . If there 
are several wells nearby, their signals are weighted accord- 
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ing to the strength of their coupling to the SET; in this case. 
Coulomb interactions between wells need also be taken into 
account. 

At T = 0, the compressibility has a (5-function peak at a 
crossing of energy levels between N and + 1 QPs in a well. 
At finite T, the peak broadens and may shift as a function of /i 
due to entropy effects. The simplest case to consider is an iso- 
lated well at the transition from iV = to = 1, or slightly 
more involved, from one to two. At higher occupation num- 
bers, we give qualitative arguments for the stability diagram. 
We examine both circular and elliptical traps, and account for 
temperature effects including broadening and excited states. 
In our model, e/4 QPs are represented as interacting charged 
particles in a magnetic field, confined to the lowest Landau 
level (LLL), with non-Abelian statistics. The interaction is 
Coulomb, supplemented by an interaction Vx (r) due to the 
exchange of Majorana fermions. 

Qualitative Picture We begin with the charging diagram 
for Abelian particles in a well, to contrast it with the non- 
Abelian case. For simplicity of presentation, we consider 
varying only the QP chemical potential, although as discussed 
below, local gating will be required to access the full charg- 
ing spectrum. The well sits in a larger quantum Hall state 
containing other distant wells, which provide a reservoir for 
QPs. At T = 0, as a function of chemical potential, a series 
of peaks in the compressibility appear, corresponding to indi- 
vidual charging events in the well. The spacing of these peaks 
defines the charging energy, U{N). As T increases, the peak 
centres evolve vertically in the charging diagram (red dashed 
lines in Fig [T}, until T reaches the minimum excitation en- 
ergy, Aat, set by the excited states within the well. Above this 
energy, the curve deviates from a straight line due to entropic 
effects. The peaks broaden linearly with T for both Abelian 
and non-Abelian QPs. 

When several non-Abelian QPs occupy a tightly confining 
well, they uniquely fuse at low energies. This produces a dis- 
tinct experimental signature - the even-odd effect. As high- 
lighted in |4|, the density dependence of the zero-temperature 
entropy produces a distinct signature in the inverse compress- 
ibility of bulk samples at low T. In local traps with dis- 
crete QP number, the difference in zero-temperature entropy 
AS between adjacent number states produces a related low- 
T signature in the charge stability diagram: the slope of the 
charge transition line in the /i — T plane is —1/AS. The 
first QP placed in the well contributes Sna = In 2/2 to the 
non-Abelian entropy (fcs = 1), or equivalently adds a \/2 
degeneracy; thus, the = — 1 transition line has slope 
—2/ In 2 in the ^ — T plane as T — > 0. A second QP fuses 
uniquely with the first QP into the 1 ot tp channel and the 
non-Abelian entropy is extinguished, AS' = — In 2/2. Thus, 
as T increases from zero, the = 1 state becomes entropi- 
cally more favourable than the N — 2 state, and the transition 
line has slope +2/ In 2 (blue solid lines in Fig[T]l. This even- 
odd effect persists as the well charges: odd numbers of parti- 
cles fuse into the non-Abelian ci-channel, while even numbers 
uniquely fuse into either the Abelian 1 or ^ channels, as long 



as T remains below the splitting between these two channels. 
If the splitting between channels is smaller than An, there 
exists an intermediate regime, in which the degeneracy is pre- 
served and the non-Abelian entropy increases by In 2 /2 with 
every additional particle, and all lines have parallel negative 
slopes. A similar effect for electrons due to spin degeneracy 
was predicted and seen in quantum dots at B = Il2l IT3\ . 

For wells far apart compared to the magnetic length, the 
rate, 5^^\\/h of Majorana exchange between them falls off ex- 
ponentially in their separation. We therefore consider charg- 
ing lines for temperatures T ^ f^weii, assumed zero for an 
isolated well. In this limit there is no fusion-channel splitting 
between wells, and each independently exhibits the even-odd 
effect. However, the charging spectra are not completely in- 
dependent due to capacitive coupling. In an experiment sen- 
sitive to multiple disorder-induced wells, the charging spectra 
of the wells appear overlaid with unknown offsets making the 
even-odd effect more difficult to observe, without first asso- 
ciating the various peaks to their respective wells. Experi- 
ments |8] [141 ESI suggest that determining such associations 
is possible. 

Equilibration Although QPs are locally trapped, the equi- 
librium model we present requires that the system explores 
the degenerate ground state manifold faster than the measure- 
ment time texp of the charge-sensing experiments. This time 
is determined by the rate of change of the gate voltage, typ- 
ically texp ^ 0.1s. We estimate the equilibration time due 
to thermal excitations as ti ^ lO^'^s ^ tgxp, meeting the re- 
quirement 1 11 1 . Moreover, the observed changes in the charge 
state of the studied well during experiments |8 1 imply that QPs 
hop freely between wells on the time scale tgxp- Assuming the 
hopping processes have a stochastic component, they will nat- 
urally lead to braiding of QPs from different wells. 

Quantitative Picture Returning to a single well in a large 
bath, we present a model for calculating the charging diagram. 
The partition function is 



Z = Y,g{N)e 



-l3{F{N,T)~piN) 



(1) 
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where /? = 1 /T, the internal free energy of A^^-particle states 
in the well is F(Ar, T), and 



9{N) = 



^2 A^odd 
1 A^ even 



(2) 



captures the non-Abelian degeneracy associated with net fu- 
sion within the well. In a well where QPs are close, such 
that all other fusion-degeneracies are split by energies larger 
than T, we take F{N,T) w F{N,Q) for T < Aat, the 
gap to excitations. In principle, however, for a wide well 
where electron-electron interactions localise the QPs further 
apart, an intermediate regime can exist in which the topolog- 
ical degeneracies are not significantly split and F{N, T) « 
F{N, 0) - T[f J In 2, where [-J denotes the integer part, for 
temperatures up to the gap A^. 

The compressibility follows from the partition function. To 
leading order near the A^ — 1 to A^ charge transition at the 
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critical chemical potential, fi^ = ^(^' 0) - -1,0), 



9{N) 0{S^-AF) 



(3) 



N 



where J^i = fi-fiS and AF = F{N,T)~F{N~l,T)~fio 
We differentiate with respect to to find the centre of the 
peak: Sfi^^^ = Tin {g{N - l)/g{N)) + AF. For a tightly 
confining circular well at low T, for which AF = 0, this 
gives (5/iinax — ±(r/2)ln2, which confirms that the charg- 
ing line slopes alternate sign as a function of the parity of N. 
In the intermediate regime, the slope is negative for all N. 
The peak height decreases with T as k^i^x ^1/ 4r, while the 
full-width-half-max (FWHM) increases with T due to num- 
ber fluctuations as FWHM - 2Tln(3 + 2\/2), roughly ten 
times as fast as the shift in position. Nevertheless, tracing 
the peak should be possible if measurements are sufficiently 
accurate. In the experimental regime of interest, the charg- 
ing energy U{N) — fi^ — fJ-Q^^ ^ T, so the peaks remain 
distinguishable. The key input to the above statistical model 
is the microcanonical low-energy spectra of fixed numbers of 
QPs in a well, which we now calculate. 

One Particle For a particle in an elliptical harmonic well. 



ay 



(4) 



where k is the spring constant, and a controls the eccen- 
tricity (a = 1 defines a circular trap), the level-spacing is 
Ai = ky/al*g^, where 1*^ = yJh/e*B is the effective mag- 
netic length for QPs in a magnetic field B. At finite T, this 
produces an internal free energy. 



F(1,T) = rin(l 



-Ai/T 



(5) 



This free energy decreases weakly with T for T < Ai, only 
significantly correcting the linear charging curve for T ^ Ai, 
as shown in figure [T] 

Two Particles As the fusion channel, /, of two orbiting 
non-Abelian anyons is conserved, the orbital dynamics may 
be treated separately in each /-sector This reduces to the 
dynamics of Abelian anyons whose statistical angle 9 de- 
pends on the fusion sector For Ising anyons, = and 
9^, ~ 7r/2 |16|. To model two such anyons in a well, each 
with charge e* = e/4, we write the Hamiltonian for a pair of 
bosons with a statistical gauge field: 



coordinates in Ui are relative to the other particle. We assume 

« 2 

that the OPs interact via a Coulomb interaction, Vj = — , 
where e = e^eo is the electric permittivity of the material. This 
approximation is valid assuming that QPs do not come within 
l*Q of each other. Vx is the direct energy splitting of the fusion 
channels due to virtual exchange of Majorana fermions. It is 
related to the fusion channel splitting discussed in |9, 10|, and 
should consist of an exponential decay and oscillations, each 
on the order of several 1*^. For circular wells, the behaviour 
of Abelian anyons has been treated previously llT7]j20l . We 
summarise key results, and include corrections due to eccen- 
tricity. 

In the symmetric gauge for harmonic traps, the centre-of- 
mass (CM) and relative (REL) coordinates decouple. In the 
CM coordinate, the statistical gauge field falls out, leaving 
a single particle projected into the LLL in a harmonic well. 
For the REL coordinate, the particle is confined to a half-plane 
with the origin removed, and remains 1 18 1. We change the 
gauge, so that = 0, giving a twisted boundary condition, 
"^JRELifjTr) = e*^^ V'iJBil^i 0). The potential landscape in 
the half -plane is defined by strong Coulomb repulsion near the 
origin together with the harmonic trap, Virap + Vj, for Vx — 0. 
The twisted periodic boundary conditions allow only angular 
momenta £ = 2n + 9 /ti, for n integer. The REL-coordinate 
wave-functions in the LLL have a basis given by \£), 



{z\t,=N^ 



(7) 



where z ~ x + iy and Ng is a normalisation constant on the 
half-plane. In this basis, we can diagonalise to find the two- 
particle spectrum. The potential has diagonal terms, as well 
as an off-diagonal term only when circular symmetry is bro- 
ken JH. 

Circular Well We assume Vx = initially, and note that 
the CM coordinate behaves just like the single particle case 
with AcM ~ Ai. The lowest energy gap A^ in the relative 
coordinate within a fusion channel / can be found by tak- 
ing differences between adjacent ^-states near the minimum, 
obtained by diagonalising the Hamiltonian. We define the pa- 
rameter To — (e*^/27refc)^^'^, the radial position of the min- 
imum of the potential. This splitting A^ oscillates with ro 
at fixed magnetic field with an amplitude that decays in the 
large-well limit, tq ^ Ig, as 



< 12Ai 
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(8) 



1 ^ 2 

^=^E(P^'-H— eM,) +Kap(rl)+ (6) 



+ V7(ri -ra) + Fi(ri -fs) 



The first term contains the electromagnetic vector potential 
Ai, corresponding to a uniform external i?-field, as well as a 
statistical gauge field (a^, aQ = (y, — x), which binds a 
flux tube of strength 0/ to each quasiparticle, and m is the ef- 
fective QP mass. We project into the LLL, taking m ^ 0. The 



The other relevant gap for the relative coordinate is the energy 
difference Ei^ — \Eq — Eq\ between lowest energy states in 
the 1 and ip channels. With Vx — 0, the splitting between fu- 
sion channels is an interchange effect, which follows from the 
allowed angular momenta in each channel; in particular, Ei^ 
behaves similarly to A^ with a maximum oscillation bounded 
by the power law ff^'^, which is approximately 20% of 

the amplitude of A^. For T < Ei^ and A^, the slope of 
the 1-2 transition in the ii — T plane is positive, exhibiting the 
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even-odd effect. Clearly, intra-channel entropy washes out the 
effect for T > A[. As 7'o varies, Ei^ will oscillate in sign, 
and can be arbitrarily small if tq is close to a zero-crossing. If 
Ei^ < T < A^, the 0-1 and 1-2 charging lines are parallel 
with negative slope —2/ In 2 ifTTl . 

Non-Abelian QPs at finite separation can exchange Majo- 
rana fermions, leading to an additional fusion channel split- 
ting. Unlike the orbital contribution, this splitting occurs even 
when QPs are localised. Using a variational method to cal- 
culate this energy splitting for particles on a sphere, it was 
found to decay exponentially on the order of several magnetic 
lengths, up to a numerical pre-factor of 0(1) [10|. Vx in 
the Hamiltonian accounts for a splitting of this form. We do 
not calculate Vx explicitly, but note that while it dominates 
the shift between fusion channel spectra in tightly confining 
wells, it oscillates and decays exponentially as the well widens 
and particle separation increases. In general, Vx increases 
Eifj,, promoting the even-odd effect over a larger T-range, and 
making a regime of parallel charging lines less likely. 

Anisotropic Well For anisotropic wells, again taking 
Vx — initially, consider the relative coordinate for two 
QPs. Starting from the circular well where QP orbits en- 
circle the origin in the half-plane, as the eccentricity a in- 
creases, the effective potential acquires a minimum on the 
X-axis, at a; = tq, and a saddle point on the y axis at 
y = ri^ja^l'^. For any given a > 1, the wavefunction be- 
comes effectively localised near the potential minimum for 



stant k — Y'^ )^ 
tum wells, B = 3.5T and d 



X{a). This is when the lowest- 



energy state near the minimum has energy lower than the sad- 
dle point potential. As a — )• 1, \{a) diverges as {a — 1)^^^^, 
confirming that for a circular well, QPs are not localised. Ec- 
centricity breaks any accidental degeneracies which arise in 
the circular potential near r^, and modifies the spectrum of 
the well. For low eccentricities, the degeneracy breaking can 
increase or decrease the orbital splitting. For large enough 
a, the QPs are trapped at opposite ends of the well, and no 
longer orbit each other, except for quantum tunneling across 
the saddle point. In a saddle point tunneling model, the or- 
bital exchange rate, R, in the large well limit is Gaussian in 
the well-size, R w Wg^ exp[-a-~^^^X{a)-^c{a) {ro/l*B)\ 
where c(a) depends weakly on a and goes to a constant of or- 
der unity as a 1 [1 1]. This expression may be obtained by 
estimating the potential as Harmonic near the minimum, and 
using a WKB type calculation of the tunneling of a particle 
near a quadratic saddle point in the LLL, as in |2r|. Increas- 
ing a also has the effect of raising the energy of the ground 
state, by increasing the harmonic frequency of the trap. 

For anisotropic wells with Vx 7^ 0, the exchange effect nat- 
urally dominates the splitting at large tq, since the exchange 
of neutral Majorana fermions decays exponentially while the 
interchange of localised charged particles in a magnetic field 
decays as a Gaussian. We recover the even-odd effect for T 
below this splitting, regardless of QP locaUsation. 

Energy Estimates A simple model producing a charge 
trap is provided by considering a point-like gate, a distance 



d above the 2DEG. A charge +|e| on this gate produces an 
effective circular harmonic trap in the plane with spring con- 
Using er ~ 13 for GaAs/AlGaAs quan- 
lOOnm, we find tq — 63nm. 
The charging energy is 1.6K, and the gap to single parti- 
cle excited states in the well is Ai « 0.24K, preserving 
the slope of —2/ In 2 throughout the accessible experimental 
range 20mK < T < 80mK. The 1-channel ground state has 
lower energy than the ^/^-channel by Ei^p w 29mK in the ab- 
sence of Vx, and the intrachannel gap Aj^ w 220mK, above 
the accessible range. As ro/l*g w 2.3, we expect the contri- 
bution of Vx to enhance the even-odd effect. Since the cal- 
culated charging energy is larger than the energy gap for the 
1/ = 5/2 plateau, it is probably impossible to observe multiple 
transitions in a single well simply by changing the voltage on 
a back gate. However, applying a voltage to a point-like gate 
on top of the sample can change the depth of a well by a large 
amount without inducing QPs in the surrounding 5/2 state. 

To further enhance the even-odd effect, all energy gaps need 
to be increased. Increasing the charge on a point-gate or re- 
ducing the setback distance d makes the confining trap tighter 
Increasing the magnetic length - by lowering B while main- 
taining the filling fraction - increases all of the relevant split- 
tings in a fixed trap geometry. 

Conclusion The detection of non-Abelian QPs through 
local charge-sensing measurements falls within realistic ex- 
perimental parameters. A sensitive compressibility measure- 
ment could extract slopes reflecting the degeneracies of the 
ground state. Additional control over confinement potentials 
will allow for even more conclusive experiments. 
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saddle point stated in the paper. In the REL coordinate for two 
particles, the potential is: 
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with ( = e*^/47re and effective magnetic length V^I*b- In the 
circular case, a = 1, the minimum is circularly symmetric at 
ro = {2(/ky/^. For a 7^ 1, in the x-y plane, (ro, 0) is still a 
minimum, and (0, r^) is a saddle point, with rg ~ r^jifii^l'^'). 
Expanding to quadratic order near the minimum and the sad- 
dle point gives: 



- ■ro,y) 



V{0,y-r,)^V{r,) + 
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where as expected the transverse components vanish in the 
circular case, and the energy levels near rg are £'„ = V{ro) + 
m*B^y'{a- l)/3{n + 1/2). We compare the ground state 
near tq to the height of the saddle point, V{rs), namely: 



SE V{rs) - Eo 



3fc 
Y 



1/3 



a 



(12) 

Solving for > gives the condition {ro/l*B)'^ > A(a) de- 
scribed in the text. To find the tunneling rate through the sad- 
dle point, we use the expression for the transmission through 
a saddle point potential Vsp{x,y) ~ —UxX^ + Uyy^ given 
in lEl, T = (1 + cxp(-7ri;))-\ where E = -5E/ei, 
ei — \/UxUy/mu]c in the large S-limit, and uoc is the cy- 
clotron frequency. For the saddle point under consideration, 
TOWc = l/2^B^ and ei = |^3a(Q; - 1)1*b^ , giving: 



E 



r2(aV3_i)_;| 



V'3a(a - 



(13) 



(\ 1/3 
jJ-J ""^j^^ for large rn/lg, 

and in this limit, T ^ cxp{ttE). To convert from transmission 
probability to a transmission rate, we take the velocity of a 
QP about its orbit as the ratio of the gradient of the potential 
to the magnetic field, and dividing by the circumference of an 
orbit near the minimum, we find the frequency of the orbit is 
~ klg'^, which, multiplied by T, gives the rate R reported in 
the main text. 



2-particle Splittings 



a = 1. 




1.0 1.1 1.2 1.3 1.4 1.5 1.' 

Hro/ts) 



1.7 1.8 1.9 



FIG. 2. log{E-i^ / Eq) calculated numerically by exact diagonalisa- 
tion, vs. log(ro//|j) for a circular well. A power law r^^ is also 
plotted as a guide to the eye. i5o = | -pr , as described in the 
paper below equation 8, with the vo/Ib dependence factored out. 
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Maximum Values 

The maximum values for the gaps A-^ and Ei^ in the two- 
particle REL spectrum are given in the text. The pre-factors 
are found by expanding the potential near the minimum at ro, 
finding the lowest-energy and first excited states in terms of 
allowed angular momenta, and expanding terms as a function 
of l*g/ro ^ 1. To find the angular momentum correspond- 
ing to the lowest energy state, we use the relation r = V^llg 
to find £o, the (possibly not allowed) angular momentum cor- 
responding to To, we first find the allowed angular momenta 
right above and below: 

= 2l{e - 6l//7r)/2j + Of/ir (14) 
£+ = 2(L(^-0//7r)/2j +1) + 0;A, (15) 

where [-J is the integer floor function. Next, convert back to 
positions corresponding to these momenta, and plug back into 
the potential to check which state has lower energy. If £g is 
the angular momentum of the ground state, then the angular 
momentum of the first excited state is £i — £g± 2. This argu- 
ment is sufficient for crudely extracting the large ro behaviour 
of the intra-channel gap A-^. A similar calculation produces 



FIG. 3. \og{Ei^ / Ec) calculated numerically by exact diagonalisa- 
tion, vs. tq/Ib for an elliptical well, with a = 1.3. A Gaussian 
decay exp (— cro) is also plotted as a guide to the eye. Ec — |^ 
is the Coulomb energy scale. 

Ei^. The decay of Ei^, is plotted in figure |2] for a circu- 
lar well; the log-log plot demonstrates a power law decay as 
Tq^. Figure [3] shows the decay of Ei^ for an elliptical well 
with a = 1.3; the quadratic decay in a log plot confirms a 
Gaussian form as expected from the saddle-point calculation. 



Relations between k , the local incompressibility, and the SET 
signal. 

Chemical Potential v. Backgate 

In this letter, we calculate k, the change in local average 
number of QPs as a function of the local chemical potential of 
the QPs. Charge stability diagrams are then drawn as a func- 
tion of the QP chemical potential. The proposed SET mea- 
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surements reveal the inverse compressibility by measuring the 
change of the local electrostatic potential in the 2DEG as a 
function of a global backgate voltage. The electrostatic poten- 
tial is then converted to chemical potential of electrons, under 
the assumption that electrochemical potential is held constant. 
A capacitance model is used to convert backgate voltage to av- 
erage electron density. The qualitative even-odd effect is in- 
dependent of definitions, but the quantitative slope predicted 
needs to be scaled to match experimental parameters. 



Relation of the SET signal to the local incompressibility 

We elaborate here on the relation between the local quasi- 
particle compressibility, discussed in this paper, and the sig- 
nal measured in an SET experiment, which is commonly de- 
scribed as measuring the local electronic incompressibility. 
Throughout the remainder of the appendix, we consider a 
sample with its top surface at z = 0. A 2DEG is found a 
distance d below the surface, and a backgate is placed a dis- 
tance D below the 2DEG, at z = -(d + Z?). The SET will be 
placed just above the sample surface, z = 0+. 

A more precise description of the SET signal is that it mea- 
sures 9$set/9Vg, where $set is the electrostatic potential 
at the SET, and Vq is the back-gate voltage. In this measure- 
ment, the electrochemical potential of the 2DEG is held fixed 
at a voltage V, by connecting it to a metallic lead. (In ref fS], 
the lead is grounded; i.e., the 2DEG is in equilibrium with a 
ground surface at infinity, and we may take V — 0.) Further- 
more, 



Moreover, if we define 7 as the spatial average of the local 
incompressibility 7(r), one finds 



5$ 



SET 



= J (frK{r 



fo)<5$(r0, 



where ^{r) is the electrostatic potential at a point r just above 
the plane of the 2DEG, inside the GaAs, fp denotes the hori- 
zontal location of the SET probe, and the precise form of the 
kernel K depends on the height of the SET probe above the 
semiconductor surface, the depth of the 2DEG, and dielectric 
constant e of the material. In general, the fluctuation in ^set 
may be interpreted as a weighted spatial average of $ within 
a distance of the order of the SET-2DEG separation, d. 

It is customary to define a local chemical potential for elec- 
trons in the 2DEG, by 

where e < is the electron charge. The change in the average 
density of electrons in the 2DEG produced by a change in the 
back-gate voltage is 

where C is an effective capacitance per unit area. Then, if 
we define the local electronic incompressibility by 7(r) = 
d^e (0 / see that 
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Sr K{f ~ fo) 7(r), 



where D is the distance to the back gate. We note that 7 will 
be finite, even when the bulk of the system sits in the quantum 
Hall plateau, due to the effects of changing quasiparticle pop- 
ulations in wells whose depths are close to a critical value. In 
practice, in the experimental geometry where D is the order 
of a micrometer, the first term will be large compared to the 
second, and C will be determined primarily by the geometric 
capacitance. 



Relation between the SET signal and the local quasiparticle 
compressibility 

As a consequence of Poisson's equation, the value of $(r) 
will be directly affected by changes in the local electron den- 
sity n(j) . In the simplest case, we consider a situation where 
there is a single chargeable potential well, surrounded by a 
region of incompressible 5/2 state, in the area sensed by the 
SET. Then changes in the local electron density result primar- 
ily from changes in N , the quasiparticle occupation number 
of the well. To a good approximation, 

5$SET , 

where 77 is the model-dependent potential at the SET from a 
single QP, and c\ is a slowly varying number accounting for 
the change in density in the rest of the sample beyond the well. 

The quasiparticle chemical potential /i employed in our pa- 
per is related in a complicated way to the local electron chem- 
ical potential /Ltg. Roughly, 

b[i = -e*(5l> 

where $ is the electrostatic potential at the position of the 
well, excluding any potential due to the presence of one or 
more quasiparticles in the well. When the expectation value of 
iV is a rapidly varying quantity, there will be a very large dif- 
ference between the variation in $ and the variation in <I>(r). 

In general, we expect that $ should depend smoothly on 
the back-gate voltage Vg- We may estimate this dependence 
by assuming that the potential well is surrounded by an in- 
compressible region of radius R, and that outside this region 
we have a continuous medium characterized by a finite in- 
compressibility 7, which we identify with the spatial average 
of 7(7^ defined above. We shall assume that D is very large 
compared to R and to the "screening length" Is = 7e/e^, but 
we should still consider different values of the ratio R/lg- In 
the case R > Ig, analysis of the resulting electrostatics prob- 
lem leads to a result 

5$ - C2(5Vg R/D 
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where C2 is a constant of order unity. In the opposite hmit. 
Is > R, we find 

<5$ - SVgls/D = SVaj/Ce"^. 

When the system is in the middle of a quantum Hall plateau, 
we expect that 7 will be large, and we might expect to be in 
the regime Is > R. 

Finally we may put these results together to find the relation 
between the SET signal and the quasiparticle compressibility 
defined in the text. In the limit where Ig > R, we have, ignor- 
ing a smooth background contribution. 



SET 



dVr, 



e' \ £7 K r]{d) 
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D 



The above model can be generalized to a situation where 
there are several wells beneath the SET tip, by choosing a 
larger radius R within which there is no continuum back- 
ground compressibility, and including explicitly the Coulomb 
interactions between quasiparticles in different wells in this 
region. In the simpler model we have replaced all wells by a 
continuum, except for the one under consideration. 



Calculation of r;(d) 

We calculate r/(c?) for a point-hke SET, located an infinites- 
imal distance above the sample surface, which is a distance d 
from the 2DEG. The distance to the backgate, D is taken to be 
much larger than d. We consider a QP of charge e* added to 
the 2DEG in a bulk sample with permittivity e, and we want to 
know the potential at the SET, when the sample sits in vacuum 



- i.e outside, permittivity is eq. We have to solve the following 
equations for an electric field E: 

eV ■ E = p, z <0 
€0^ ■ E = 0, z>0 
V X = 0, everywhere, 

with boundary conditions at the sample boundary (i.e. z = 
0) of continuous E'-fields in the x and y directions, and 
lim cqE^ = lim eE^- Place an image charge q' at z = 

d above the surface, and then using cylindrical coordinates 
(r, (j), z), the potential at any point inside the sample is: 
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7re 



R~i 



R2 



z < 0, 



(16) 



where i?i = ^/r^ + (d + z)^, R2 = v^r^ + (d-z)^. For 
the region z > 0, which is where the SET is, there are no 
charges, and the potential must therefore be a solution to 
Laplace's equation without singularities. The simplest solu- 
tion is the potential from an effective charge q located at the 
site of the QP e*, giving a potential: 



1 q 



z>0. 



47reo Ri 

The solutions [16] and [17] can be matched at z 
satisfy the boundary conditions, giving e* — q' 
q')/e ~ q/^Q- This implies q = (2eo/e + eg) e*. We thus find 
the potential at the SET due to the QP to be: 



(17) 

and must 
q and (e* + 



77(d) = 



1 



47r e + eq d 

A very similar calculation can be found in section 4.4 in |22|. 



